Computing Effective Pair Potentials from Pair Distribution Data

c l a s s M C A b s t r a c t O b j e c t { p u b l i c : / / D o a l o c a l m o n t e c a r l o s t e p . / / R e t u r n t r u e i f t h e s t e p w a s a c c e p t e d , f a l s e o t h e r w i s e . v i r t u a l b o o l s t e p ( ) = 0 ; } ; / / A b s t r a c t b a s e c l a s s f o r M C r u n s . c l a s s M C A b s t r a c t R u n { t y p e d e f v e c t o r < M C A b s t r a c t O b j e c t * > C o n f i g u r a t i o n ; p r o t e c t e d : C o n f i g u r a t i o n c o n f i g u r a t i o n ; p u b l i c : / / c o n s t r u c t o r M C A b s t r a c t R u n ( ) : c o n f i g u r a t i o n ( ) { } / / d e s t r u c t o r v i r t u a l ~ M C A b s t r a c t R u n ( ) { } / / a d d o n e o b j e c t t o t h e c o n f i g u r a t i o n v i r t u a l v o i d a d d _ t o _ r u n ( M C A b s t r a c t O b j e c t & o b j e c t ) { c o n f i g u r a t i o n . p u s h _ b a c k ( & o b j e c t ) ; } / / r u n t h e s i m u l a t i o n v i r t u a l v o i d r u n ( ) { w h i l e ( ! e n d _ c o n d i t i o n ( ) ) { d o _ m c _ s t e p ( ) ; e v a l ( ) ; } } / / d o o n e M C − s t e p v i r t u a l v o i d d o _ m c _ s t e p ( ) { u n s i g n e d i n t l o c a l _ s t e p ; f o r ( l o c a l _ s t e p = 0 ; l o c a l _ s t e p < c o n f i g u r a t i o n . s i z e ( ) ; l o c a l _ s t e p + + ) / / d o a l o c a l s t e p i n a r a n d o m p a r t i c l e c o n f i g u r a t i o n [ m c _ r a n d ( c o n f i g u r a t i o n . s i z e ( ) ) ] − > s t e p ( ) ; } / / d o t h e e v a l u a t i o n v i r t u a l v o i d e v a l ( ) = 0 ; v i r t u a l b o o l e n d _ c o n d i t i o n ( ) = 0 ; / / * * M e t h o d s t o u s e d u r i n g e v a l u a t i o n v i r t u a l c o n s t C o n f i g u r a t i o n & g e t _ c o n f i g u r a t i o n ( ) { r e t u r n c o n f i g u r a t i o n ; } } ; # e n d i f M C A bs tr ac tR un .H # i f n d e f M C R E L A T I O N C O U N T E R _ H # d e f i n e M C R E L A T I O N C O U N T E R _ H / / T h i s c l a s s p r o v i d e s a s i m p l e a c c e p t a n c e − r a t e c o u n t e r t o b e u s e d / / w i t h M C O b j e c t s . S i m p l y c a l l t h e c o u n t _ a c c e p t e d ( ) / c o u n t _ r e j e c t e d ( ) / / f r o m t h e a c c e p t ( ) / r e j e c t ( ) − m e t h o d o f t h e o b j e c t a n d g e t ( ) r e t u r n s / / t h e a c c e p t a n c e r a t e s i n c e t h e l a s t c a l l o f r e s e t ( ) . c l a s s M C R e l a t i o n C o u n t e r { p r o t e c t e d : i n t c 1 , c 2 ; p u b l i c : M C R e l a t i o n C o u n t e r ( ) { r e s e t ( ) ; } v i r t u a l ~ M C R e l a t i o n C o u n t e r ( ) { } v o i d r e s e t ( ) { c 1 = 0 ; c 2 = 0 ; } o p e r a t o r d o u b l e ( ) { r e t u r n g e t ( ) ; } d o u b l e g e t ( ) { i f ( c 1 + c 2 = = 0 ) r e t u r n 0 . 0 ; r e t u r n ( d o u b l e ) c 1 / ( d o u b l e ) ( c 1 + c 2 ) ; } v o i d i n c 1 ( ) { c 1 + + ; } v o i d i n c 2 ( ) { c 2 + + ; } } ; # e n d i f M C R el at io nC ou nt er .H 2/ 22 ba se /M C A bs tra ct R un .H , u til s/ M C R el at io nC ou nt er .H